Ultrafast dephasing of light in strongly scattering GaP nanowires. 
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We demonstrate ultrafast dephasing in the random transport of light through a layer consisting 
of strongly scattering GaP nanowires. Dephasing results in a nonlinear intensity modulation of 
individual pseudomodes which is 100 times larger than that of bulk GaP. Different contributions to 
the nonlinear response are separated using total transmission, white-light frequency correlation, and 
statistical pseudomode analysis. A dephasing time of 1.2 ±0.2 ps is found. Quantitative agreement is 
obtained with numerical model calculations which include photoinduced absorption and deformation 
of individual scatterers. Nonlinear dephasing of photonic eigenmodes opens up avenues for ultrafast 
control of random lasers, nanophotonic switches, and photon localization. 

PACS numbers: 42.25Bs, 78.47J-, 78.67Uh 



The appearance of an interference pattern after trans- 
port of coherent light through a multiple scattering 
medium is the result of coherent summation of thousands 
of light paths with random phases. In recent years, new 
methods have been developed to exploit the coherent as- 
pects of diffuse light transport for imaging [l], Q . Clas- 
sical waves, like light or sound, provide unique opportu- 
nities to study mesoscopic wave transport, such as local- 
ization, in the absence of many-body interactions usually 
found for electrons 0, Q • Coherent light scattering has 
recently lead to exciting new directions of research, such 
as coupling of localized states with quantum emitters Q , 
transverse localization, and combined effects of localiza- 
tion and nonlinearity [6|. 

Although the lack of many-body interactions makes 
phase coherence generally more robust for optical waves 
than for electrons, several dephasing processes can in- 
fluence the transport of light, such as magnetic fields, 
scattering from moving particles, and polarization effects 
The time required for electrons or photons to loose 
their coherence due to inelastic collisions is known in the 
theory of electronic conductance as the phase breaking 
time [10]. While for electrons phase breaking (deco- 
herence) processes have only been measured indirectly 
through changes in the conductivity with temperature, 
dephasing can be observed dir ectly for optical waves by 
the changes in speckle pattern [11[, which is used in dif- 
fusing wave spectroscopy [l2| . 

In this Letter, we demonstrate a new regime of re- 
versible dephasing in a random medium on ultrafast time 
scales. The ultrafast dephasing takes place a million 
times faster than other dephasing mechanisms such as 
Brownian movement, opening up new applications in ul- 
trafast control of random lasers, quantum emitters, and 
localization phenomena. The material under study con- 
sists of a layer of semiconductor nanowires, which are 
important for novel applications in optoelectronics and 
photovoltaics 13-[l5|. We have recently demonstrated 



that these semiconductor nanowires form one of the most 



efficient light trapping layers available to date [16j. Here, 
we combine these favorable scattering properties with the 
intrinsic nonlinearity of the semiconductor host material 
to achieve active control on ultrafast time scales. 

Gallium phosphide (GaP) nanowires were grown us- 
ing metallo-organic vapour phase epitaxy as described in 
Ref. 16]. The length of the nanowires was controlled 



by the VLS growth time, while the thickness and vol- 
ume fraction was independently tuned by switching to a 
lateral growth regime where material was deposited ho- 
mogeneously over the nanowire sidewalls. In the experi- 
ments of this work, we make use of two nanowire samples, 
respectively with mean free path £ = 0.16 ± 0.02 /im and 
thickness L = 1.5 /im (sample 10 of Ref. [TBI], here sam- 
ple 1), and L = 4.0 /im, £ = 0.2 ± 0.02 /im (sample 11 
of Ref. [16|, here sample 2), at a wavelength of 632 nm. 
The variation of £ with wavelength was found to be small 
over the spectral range of our experiments [16 ]. 

We investigated the photoinduced nonlinear response 
using a broadband white-light supercontinuum probe 
generated by focusing of 800-nm laser pulses from a re- 
generative amplifier (Coherent RegA, 250 kHz, 200±10 fs 
pulse duration) in a quartz window. The material is 
pumped using the second harmonic at 400 nm wave- 
length. Both the pump and probe are focused onto the 
sample using a lens with a 6-cm focal length, resulting 
in a focus of around 5 /im in diameter and a pump fiu- 
ence of 8 mJ/cm 2 . Total transmission measurements are 
shown in Fig. [TJa) for the nanowires as well as for a GaP 
substrate (thin line, blue). 

In absence of absorption, the transmission through a 
scattering layer follows Ohm's law T ~ £/L. The total 
transmission can be derived including effects of internal 
reflection 17], which for small absorption (L a > L > £) 
reads 
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FIG. 1: (color online) (a) Total transmission spectra for 
nanowire layers 1 and 2, and for a GaP substrate (thin line, 
blue). Inset: scanning electron microscopy images, scale bars 
1 /mi. (b) Pump-probe total transmission AT/T for the sam- 
ples of (a) and for a 51-nm diameter nanowire layer (dash, 
sample 6 of Ref. [3). (c) Absorption length a -1 obtained 
from (b) using Eq. (pQ). 

where L denotes the slab thickness, z e ^ are extrapolation 
lengths which represent the effects of internal reflection 
at the front and rear interfaces of the diffuse slab. The 
diffuse absorption length is defined as L a = (I /3a) 1 / 2 , 
with a -1 the bulk absorption length. 

Pump-probe total transmission measurements are 
shown in Fig. QJb). All samples show a reduction of 
the total transmission directly following excitation by 
the pump laser. Independent experiments were per- 
formed verifying that mesoscopic contributions to the 
total transmission, i.e. fluctuations caused by a finite 
number of transmission channels [18|, were unimportant 
under the present experimental conditions. Therefore the 
signals in Fig. DJb) are attributed to changes in the ab- 
sorption coefficient of the materials. Transient absorp- 
tion in semiconductor nanostructures results from free 
carrier generation [l9|. A faster relaxation time is ob- 
served for nanowires compared to the bulk semiconduc- 
tor, which can be attributed to trapping of photogener- 
ated carriers at surface states [20j . The transient rela- 
tive absorption increases from 3 x 10 -3 for the GaP sub- 
strate to 5.2 x 10 -2 for nanowire sample 2. We attribute 
this strong increase in absorption by light trapping in 
the nanowire layer [2l[. Light trapping can be assessed 
from the time the light resides in the micrometer-thin 



photoexcited region, where it is susceptible to nonlinear 
absorption. For the GaP substrate, the coherent beam 
traverses this region in around 10 fs, while, for nanowire 
sample 2, the light diffusion time amounts to approxi- 
mately L 2 /tt 2 D c± 0.15 ps. We note that optically thin 
nanowires [dashed line, green in Fig. OJb)] did not show 
such a marked increase of the nonlinear absorption, elim- 
inating other possible nanowire surface effects. 

Equation [T] can be inverted to yield values for the 
absorption length a -1 for a pump-probe signal AT/T, 
as shown in Fig. DJc). A minimum absorption length 
a -1 of 49.9 ± 4 jam is found, corresponding to a change 
of the imaginary refractive index Aft = coa/47ru of 
(1.0 ± 0.2) x 10 -3 . Using a free-carrier absorption cross 
section of GaP of 1.1 x 10 -18 cm -2 [22[, our value for 
Aft amounts to a free-carrier density of 2 x 10 20 cm -3 , 
which is a realistic density for the strong optical pumping 
regime under study. 

In the following we focus on the results for nanowire 
sample 2. Similar results are obtained for sample 1 (23| . 
To assess the effects of nonlinear pumping on the dy- 
namical transport parameters of the random medium, we 
measured the frequency correlations of light transmitted 
through the nanowire slab. For this purpose, a cone of 
light of around 2° angular width was collected and ana- 
lyzed using a fiber spectrometer. Broadband linear po- 
larization filters were used both for the incident and for 
the detected light. Spectra were measured as a function 
of pump-probe delay time, yielding time-resolved spec- 
tral maps such as shown in Fig. Efa). The transmission 
spectra consist of large fluctuations corresponding to fre- 
quency speckle. Within 1 ps following excitation with 
a pump laser, changes occur in both the spectral posi- 
tion and amplitude of the fluctuations, as is illustrated 
in Fig. [2fb) where we compare two spectra taken respec- 
tively 3 ps before and 1 ps after arrival of the pump pulse. 
To quantify these variations in the speckle pattern, we 
calculated the time-correlation function C a h(y \ t, t') from 
the time-resolved spectra, using the conventional defini- 
tion of the intensity- intensity correlation [24|. The spec- 
tral cross-correlation was calculated for each spectrum 
at time t with the spectrum at t' = — 3 ps, as shown in 
Fig. [2^c), where spectra taken at 20 different locations 
on the sample were used to obtain an ensemble average. 

For times directly after the arrival of the pump pulse, 
the time-correlation shows a reduction in amplitude, ac- 
companied by a spectral shift, denoted hereafter by 5vo, 
of 10.2 ±1.0 cm -1 . Cross sections of Fig. [21(c) are shown 
in Fig.[3](a) at -3 ps (black dots) and 1 ps (red diamonds). 
The spectral autocorrelation can be calculated using [24[ 

C ab (Av, Ta ) ~ N[\ V L\ 2 /\smh( V L)\ 2 } , ( 2 ) 

with rj = D~^(t~ x +27riA^) 1 / 2 , where D is the diffusion 
constant of the light, and r a = (vEOt)~ l the absorption 
time, ve = 3D/£ being the energy velocity. The nor- 
malization constant N is used to fit Eq. ([2|) to our data 
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FIG. 2: (color online) (a) Map of time-resolved frequency 
speckle for sample 2, with (b) spectra taken at — 3 ps (black) 
and 1 ps (red) pump-probe time delay. Blue arrows indicate 
peaks selected for pseudomode analysis, (c) Map of experi- 
mental correlation C a h(y\ t, t'). 



given the reduced experimental speckle contrast. Equa- 
tion (j2j) gives good agreement for the autocorrelation at 
-3 ps [line, black in Fig. 02a)] for D = 15 ± 2 m 2 /s. 

The reduction of the correlation between spectra at 
different times allows for an estimate of the dephasing 
time. As shown in Ref. 8, dephasing processes affect 
the average probability of light transport through an ex- 
ponential decay exp(— t/r 7 ) of the correlation between 
two amplitudes, where r 7 is the dephasing time. The 
effects of dephasing and average refractive index change 
[iBj can thus be introduced into the correlation function 
by replacing r] in Eq. (J2j) by the new 
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We note that the normalization factor N only depends on 
the variance of the two individual configurations, which 
is not affected by the dephasing or the effective index 
change. We fitted the reduction in correlation ampli- 
tude by the dephasing model of Eq. (j3j). Good quan- 
titative agreement is found [thick line, red in Fig. [3^a)] 
for a dephasing time r 7 = 1.2 ± 0.2 ps and a change in 
the average refractive index of the scattering medium of 
An e /n e = (-5.8 ±2) x 10" 4 . 

We considered, next to the loss of correlation, also 
the effect of ultrafast excitation on individual pseudo- 
modes. Pseudomodes are defined as the transmission 
modes of the system, which take on the form of spec- 
tral resonances, with a spectral width determined by the 
leakage of the state to free space. The spectrum of pseu- 
domodes plays an important role in theories of localiza- 
tion [Isl, [I?} and random lasers 28]. We have traced the 
time evolution of a subset of pseudomodes selected from 
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FIG. 3: (color online) (a) Cross sections of Fig. [2c) at delay 
times of —3 ps (dots, black), and 1 ps (diamonds, red). Lines 
represent fits using Eq. ©, without (black) and with (red 
line) dephasing model with r 7 = 1.2 ± 0.2 ps. (b) Finite- 
element calculations [23[ of autocorrelation (line, black) and 
cross-correlations for a refractive index change AriGaP = — 2 x 
10 -3 + 10 _3 i (dashed line, blue), and an additional 10 -3 size 
increase of individual scatterers (thick line, red). 



our experimental spectra using an algorithm [29j , includ- 
ing modes such as those indicated by the blue arrows in 
Fig. Efb). The position and amplitude of each pseudo- 
mode was tracked as function of the pump-probe delay 
time, yielding a set of well-defined intensity changes and 
spectral shifts shown in Fig. 2Ja,b). Each time delay t 
represents a histogram, such as is shown in the right-hand 
panels of Fig. HJa,b) for t = 1 ps. A wide distribution of 
spectral and intensity variations is found with both de- 
creasing and increasing values, which extend beyond the 
variation of the time-correlation in Fig. [2fc). A fraction 
of modes in Fig. H^b) shows a negative frequency shift, 
which cannot be explained by an average refractive index 
change but which indicates a more chaotic change of the 
pseudomode spectrum. 

We have investigated possible dephasing processes us- 
ing a numerical finite-element (COMSOL) model of light 
propagation in a two-dimensional slab consisting of a 
random arrangement of dielectric cylinders (23[. The 
frequency cross-correlation was calculated between two 
configurations corresponding to the situations before and 
after excitation. Consistent with earlier observations 
[iBj . no decorrelation was found for only a small refrac- 
tive change of around 10 -3 . However, our model cal- 
culations identified two other possible dephasing mecha- 
nisms: absorption and particle deformation. While ab- 
sorption [dashed line, green in Fig. [3fb) for AncaP = 
—2 x 10 -3 + 10 — 3 i] does not change the phase of individ- 
ual light paths, it does lead to a redistribution of phasor 
weights in the total speckle field by suppression of contri- 
butions from long paths. This interpretation is only par- 
tially correct in time-resolved experiments, where paths 
are separated in time and only the finite duration of 
the probe pulse leads to mixing of paths with different 
amounts of absorption. We point out that inhomoge- 




I (J I ■ 1 ■ L 

t 5 10 40 20 

Time (ps) # pseudomodes 



FIG. 4: (color online) Results of pseudomode analysis of 202 
modes, yielding nonlinear intensity change I/I-3ps (a) and 
spectral shift 5iso (b). (Red dashed lines) statistical average 
of the distributions. Arrows indicate time location of the his- 
tograms shown in right-hand panels. Lines indicate averages 
(solid) and standard deviations (dash) of the distributions. 



neous absorption may also lead to variations between 
different paths of equal length and thus to dephasing. 

We explored another possible dephasing mechanism 
caused by stress- induced deformations. It was found that 
a 10 ~ 3 increase of the scatterer size is sufficient to achieve 
quantitative agreement with experiment [thick line, red 
in Fig.[3fb)]. Such deformations are likely to occur under 
the typical conditions of high power pulsed laser excita- 
tion of solids [30| , and may be enhanced by the particular 
morphology of the nanowires [311 • The effect of deforma- 
tions is important as it represents the contribution of 
changes in the scattering matrix of individual scatterers. 
No dephasing was observed when randomly displacing 
each scatterer over 0.1 nm. 

Our analysis of pseudomode distributions and spec- 
tral correlations reveals that ultrafast dephasing enables 
reversible switching between random configurations in- 
volving thousands of scattering events. Although the in- 
tensity and frequency changes appear chaotic, they are 
highly reproducible over millions of pump-probe cycles. 
The observed pseudomode intensity modulations of up 
to 50% are large for the field of nanophotonics, and com- 
parable to e.g. nonlinear switching of three-dimensional 
photonic crystals [32j . The similar decorrelation values 
for the two nanowire layers with different thicknesses [23| 
indicates that dephasing is limited by the pump absorp- 
tion length, rather than by sample thickness; further op- 
timization of dephasing using e.g. multiphoton excitation 



[32| , will be of importance for exploitation of these effects 
in applications. Deterministic control of pseudomodes 
may be combined with active gain media to produce novel 
types of ultrafast random lasers where the lasing modes 
can be manipulated using an external control pulse. The 
dephasing time of ~ 1.2 ps is of the order of the transport 
time through the layer, which paves the way for nonadi- 
abatical control over light paths. Of particular interest 
here will be the possibility to break the phase coherence 
of correlated paths involved in mesoscopic transport, and 
ultimately of time-reversed light paths involved in photon 
localization. 
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